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ABSTRACT 

In this paper we compare outcomes of some extended phantom-like cos- 
mologies with each other and also with ACDM and ADGP. We focus on the 
variation of the luminosity distances, the age of the universe and the deceler- 
ation parameter versus the redshift in these scenarios. In a dynamical system 
approach, we show that the accelerating phase of the universe in the f{R)-DGP 
scenario is stable if one consider the curvature fluid as a phantom scalar field 
in the equivalent scalar-tensor theory, otherwise it is a transient and unstable 
phenomenon. Up to the parameters values adopted in this paper, the extended 
F{R, (^)-DGP scenario is closer to the ACDM scenario than other proposed 
models. All of these scenarios explain the late-time cosmic speed-up in their 
normal DGP branches, but the redshift at which transition to the accelerating 
phase occurs are different: while the ADGP model transits to the accelerating 
phase much earher, the F{R, (p)-DGP model transits to this phase much later 
than other scenarios. Also, within the parameter spaces adopted in this paper, 
the age of the universe in the /(i?)-DGP model is larger than ACDM, but this 
age in F{G, (^)-DGP is smaller than ACDM. 

Subject headings: Braneworld Cosmology, Phantom Mimicry, Dynamical System 
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Introduction 



One of the most remarkable discoveries of the past decade is that the universe 
is currently in an accelerating phase, which means that expanding velocity of the 
universe is increasing. This phenomenon is supported by data related to the luminosit y 



measurements of high red shi ft super novae (jPerlmutter 



1999 



Riess 



1998 



Astier et al.l 



2006 



Wood-Vasey et al. 



20071 ) ■ measurement of degree-sca 



microwave background fCM B) (iMiller et al. 



large scale structure (LSS) (ICoUess et al. 



Springel et al. 



1999 



2001 



Hananv 



e anisotropies in the cos mic 



2000 



Tegmark et al 



Spergel et al. 



2004; 



Cole et al. 



2003) and 



2005 



20061 ). The rigorous treatment of this phenomenon can be provided in 
the framework of general relativity. In the expression of general relativity, late time 
acceleration can be explained either by an exotic fluid with large negative pressure 
that is dubbed as dark energy in literature, or by modifying the gravity itself which is 
dubbed as dark geometry or dark gravity pro posal. The first and simplest candidate o f 



dark energy is the co smological constant, A ( iSahni fc Wang 



2000 



Padmanabhan 



2003 



Copeland et al. 



20061 ) . But there are theoretical problems associated with it, such as its 



unusual small nume rical value ( the fine tuning problem), no dynamic al behavior and 



(Weinberg 


1989; 


Carroll 


2001: 



Caldwell et al. 



19991 ). These problems 



have forced cosmologists to introduce alternatives in which dark energy evolves during the 
universe evolution. Scalar field models with their specific features provide an interesting 
alternative for cosmological constant and can reduce the fine tuning and coincidence 
problems. I n this respect, several candidate rnodels have been proposed: "quintessence" 



scalar field ( 



Ratra fc Peebles 



1988 



Saini et al 



2000 



Brax &: Martin 



2000 



Barreiro et al 



2000l : ISahni fc Wang 



fields ( 



Caldwell 



2002 



Movassari fc Se 



Amendola et al. 



are 



2000; Sahni et al. 2002 



Tsuiikawa fc Sami 



2003 



Sami fc Padmanabhan! 120031). phantom 



2004; 



20091) and chaplygin gas (iKamenshchik et al. 



Roos 



2007 



Caldwell fc Lindei 



2005 



Bouhmadi-Lopez fc Lazkoz 



Cai et al. 



2001 



2007 



Dev et al. 



Zhang et al 



2010; 



2003 



2006al : 
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Zhane & Z 


lU 


2006b 


: Bertolami e 


t al. 


2004 




Biesiac 


a et al. 


2005 


2008; 


Roos 


2008a 


b; 


Zhane et al. 


2009 


: Setare 


2009 


) are among 



As an alternative for dark energy, modification of gravity can be accounted for the 
late time acceleration. Among the most popular modified gravity scenarios which may 



successfully describe the cosmic speed-up 



2010 



Noiiri &: Odintsov 



2003 



2004 



2005 



is f ( R) gravity (Capozziello 



2006 



Noiiri et al 



2007 



2003; 



Sotiriou 



Noiiri &: Odintsov 



2008 



Bamba et al 



2008 



2008 



Carrol 



Saavedra fc Vasquez 



et al. 



2009 



2004 



Amendola et al 



Setare 



2007 



Nozari 



2008a 



Atazadeh et al 



20081 ). Modified gravity also can be achieved 



by extra-dimensional theory in which the observable universe is a 4-dimensional brane 
embedded in a five-dimensional bulk. Dvali- Gabadadze-Porrati (DGP) model is one of 
the extra-dimensional models that can describe late-time acceleration of the universe in 
-a ccelerating branch due to 



its se. 



2000a 



b 



Dvali &: Gabadadze 



2001 



eakage of gravity to the extra dimension 



Dvali et al. 



2002 



Deffayet 



2001 



Lue 



20061) 



(iDvali et al. 



Recent observations constrai n the equation o: 



be ujx ~ —1 and even ux < —1 (IMelchiorri et al 



state parameter of the dark energy to 



2003 



Riess et al. 



2004 



Komatsu et al. 



20091 ). One of the candidates for dark energy of this kind is the phantom scalar 
field. This component has the capability to create the mentioned acceleration and its 
behavior is extremely fitted to observations. But it suffers from problems; it violates 
the null energy condition and its energy density increases with expansion of the universe 
which is an unphysical behavior. Also it causes the quantum vacuum instabilities. So, 
cosmologists have tried to realize a kind of ph antom-like behavior in the cosmolqgical 



mode 



2000 



s without introduction of phantom fields (ISahni fc Shtanov 



2003 



Sahni 



Lue fc Starkman 



2004 : 



Shtanov 



20041 ). With phantom mimicry ( the phantom-like behavior), 
effective equation of state parameter of dark component remains less than —1 and effective 
energy density of the unverse increases with cosmic expansion. In the framework of 
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general relativity these two expressions are equivalent. However, in the effective picture 
of cosmological models, in order to satisfy phantom-like behavior, the mentioned two 
expressions must be investigated separately The phantom mimicry discussed in this study 
has a geometric origin. In this paper, we discuss briefly DGP-inspired theoretical models 
that realize phantom-like behavior. We present cosmological dynamics in each proposed 
model and then we compare these models with ACDM and also with each other through 
investigation of their expansion histories. Within a dynamical system ( phase space) 
approach, we show that the accelerating phase of the universe in the /(it!)-DGP scenario 
is stable if one consider the curvature fluid as a phantom scalar field in the equivalent 
scalar-tensor theory, otherwise it is a transient and unstable phenomenon. As another 
important probe, we study the age of the universe in each model. We show that some of 
these models account for a transient accelerated phase. 



As we have pointed out in the introduction, the observational evidences show that the 
equation of state parameter of dark energy in the universe can be less than —1. Many 
attempts have been made to find dark energy models that allow the so-called phantom 
dark energy: a dark energy component with large negative pressure and negative kinetic 
energy with equation of state parameter less than —1. A phantom field is described by the 
following action 



The kinetic energy term of the phantom field in the corresponding lagrangian enters with 
opposite sign in contrast to the ordinary matter and this distinguishes phantom field from 
ordinary (canonical) fields. But this distinctive property of the phantom field causes a 
series of quantum vacuum instabilities. Also, due to characteristic p < —p feature of this 



2. The phantom-like evolution 




(1) 
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kind of matter, it violates the null energy condition. On the other hand, dark energy with 
a; < — 1 (phantom energy) is beset with a host of undesirable properties which makes this 
model of dark energy to be unattractive. So, cosmologists attempted to find a way which 
removes these problems. They are forced to consider models that realize phantom-like 
behavior without introducing any phantom matter. In the language of general relativity, 
phantom-like behavior means that the equation of state parameter of dark energy is less 
than —1 or its energy density increases with expansion of the universe. In the cosmological 
models that there is no phantom matter, two mentioned expressions for the phantom-like 
behavior are not equivalent and should be investigated separately. 



3. Phantom mimicry: extended models 



3.1. ACDM vs ADGP 



A cosmological model which has the capability to realize the phantom-like behavior 
without introducing any phantom matter, is the ADGP model. Firstly, we introduce briefiy 
the ADGP model. The DGP model, as a modified theory of gravity that modifies the 



geometr ic sector of the Einstein fie 



follows ( iDvali et al 



2000a; 



Deffayet 



d equations, ha s a modified Friedmann equation as 



2001 



Lue 



20061) 



SttG 



pit) 



(2) 



where H{t) = ^ is the Hubble parameter. This model, in its self-accelerating branch 
(corresponding to the negative sign of equation (2)), explains late-time accelerated 
expansion of the universe. But, the normal branch (positive sign), has no self-accelerating 
behavior. Nevertheless, the self-accelerating branch has beset with ghost problem. The most 
trivial extended-DGP model which can create the accelerated phase in the normal branch 
is the ADGP model. In this model a cosmological constant which plays the role of dark 



- 7- 



energy, lies on the brane. This scenario based on the so-called dynamical screening of the 
brane cosmo logical constant, has the capacity to realize a phantom-like behavior without 
introducing any phantom matter neither on the brane nor in the bulk. In this case, brane 
cosmological constant as a dark ener g y componen t is screened due to curvature modifica tion 



at late-time ( ISahni &: Shtanov 



2003 



Sahni 



2004; 



Shtanov 



2000 



Lue Sz Starkman 



20M) 



The action of this model is 



S = ^Ml j d^x^^n + ^mj j d^x^^R + j d^x^/^Cr, 



(3) 



where M5 is the five-dimensional fundamental scale. The first term in S is the Einstein- 
Hilbert action in five dimensions for a five- dimensional metric qab with Ricci scalar TZ. The 
metric q^p is the induced (four-dimensional) metric on the brane, and q is its determinant. 
The second integral contains an induced curvature effect that appears in DGP setup due 
to quantum corrections via interaction of the bulk graviton with matter on the brane. To 
apply correct boundary conditions, the Gibbons- Hawking term containing the extrinsic 
curvature of the brane should be added to the brane part of the action. The crossover 
distance, that gravity in the scales larger than it appears to be 5- dimensional, is defined as 
follows 



2M| 



(4) 



The cosmology on t he brane for a spa t ially-fiat unive rse 



Friedmann equation (jPvali et al 



2000a 



Deffayet 



2001 



Lue 



' oUows Deffayet's modified 



20061 ) 



H 



8nG 



{Pm + Pa) ■ 



(5) 



The negative and positive signs in this equation represent the self-accelerating and normal 
branches respectively. In which follows, we focus on the normal branch, because this branch 
has the key property that brane is extrinsically curved so that shortcuts through the 
bulk allow gravity to screen the effect of the brane energy-momentum contents at Hubble 
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parameter H ^ r^^ and t 



Sahni 



2004 : 



Shtanov 



l ere is no ghost instabi. 



2OOOI : 



Lue fc Starkman 



20041 ) 



ity in this branch (jSahni fc Shtanovl 12003 



^' + - = ^(Pm + PAJ. (6) 
Comparing this equation and the standard Friedmann equation in a spatially fiat universe 

= ^(Pm + PDS), (7) 



one can find an effective da rk energy component as follow (IMelchiorri et al. 



2003 



Riess et al. 



2004 



Komatsu et al. 



20091 ) 



^Pli = A - -■ (8) 

Since in this case H{t) is a decreasing function of the cosmic time, the effective dark energy 
component increases with time, and therefore we realize a phantom-like behavior without 
introducing any phantom matter that violets the null energy condition and suffers from 
several theoretical problems. We can define a Wejf by using pejj = Po(l + z)^*^^^"^'^^^-' and 



comprehend that in a red-shift such as z* w. 
energy is dominated ) , Wpf f is less than 



lich H(z*) < Ar, 



Alam fc Sahni 



2006 



Lazkoz fc Majerotto 



f|Lazkoz et al. 



during the epoch which dark 



2006 



Maartens fc Majerotto 



2006 



20071 ). In other words, in the mentioned red-shift. 



two expressions of the phantom-like behavior are equivalently applicable. Also, existence of 
^ term causes the cosmological constant A to be screened, and to be appeared less than its 
actual value. This case is called the gravitational screening effec t and as we have pointe d 



out, one can 


realize phantom-like behavior 


jased 


Sahni 


2004 




Shtanov 


2OO0I: 


Lue & Starkman 


2004) 



jased on this effect (ISahni fc Shtanov 



2003 



20041 ). The ADGP model is an alternative 



for AC DM ( the cosmological constant A plus the cold dark matter ). It has been shown 
recently that ADGP in some respects gives even b etter fit to observation s than other dark 



energy scenarios, albeit with one more parameter (jSoUerman et al. 



20091 ). To investigate 



expansion history of ADGP and comparing it with constant-w dark energy models (such 
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as the ACDM), we study variation of luminosity distances versus the redshift in these 
scenarios. By rewriting the Friedmann equation (6) in the more phenomenological form, we 
find 

^ = ^^rrXl + Zf + ^]A + - (9) 
-no 

Where 

_ SttG n _ ^ n _ 1 







Also, by using this dimensionless equation, one can imply a constrain on the parameters of 
the model as follows 

nm + nK = l + 2^/n^^. (11) 



This constraint means that a 
{VLm., ^a) plane ( 



iat ADGP model mimics a closed ACDM mode 



Laz 



Lazkoz fc Majerotto 



coz et al. 



2006 



Maartens fc Majerotto 



2006 



Alam fc Sahni 



in the 



2006 



with w 



20071 ). Friedmann equation in the ACDM model ( constant-w model 
1 ) is phenomenologically written as follows 
H(z) 







;i-n„)(i + ^)3(i+-). 



12) 



By considering luminosity-distance in the standard form in spatially flat cosmology 

Ho 



dr(z) 



l + z) 



Hiz) 



-dz , 



(13) 



we can compare dilz) calculated in ADGP and ACDM frameworks. This comparison is 
shown in figure 3. Up to parameters values used to plot this figure, the luminosity distance 
calculated in the ADGP is larger than that of the ACDM for a given redshift. 



4. /(i?)-DGP scenario 



As an extension of the previous subsection, here we consider po ssible modifica tion 



of the induced gravity on the brane in the spirit of f{R) theories (INozari 



2008a 
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Atazadeh et a] 


2008; 


Saavedra & Vasauez 


2009; 


Nozari 


2009a 


b 


d, 


Atazadeh & Sepanei 


2006, 


2007. 


2009 


). We assume that induced 



Bouhmadi-Lopez 



2009 



be modified by a general f{R) term. It has been sho wn that AD .f(R) theories can follow 



close ly the expansion history of the ACDM universe (IHu fc Sawicki 



2007 



Martinelli et al. 



20091 ). Here we study an extension of f{R) theories to a DGP braneworld setup. By 
focusing on the luminosity distance-redshift relation, we compare expansion history of this 
type of model with other alternative scenarios. The action of this model can be written as 



follows 



where by definition 



S 



Ml 



(14) 



(15) 



By calculating the bulk-brane Einstein's equations and using a s patially flat FR W line 



element, the following modifled Fried r nann equation is obtai ned (INozari 



2009a b 



Atazadeh fc Sepangi 



2006, 



2007, 



20091: 



Bouhmadi-Lopez 



20091) 



3 rr_ 



(16) 



(17) 



where 

is energy density corresponding to the curvature part of the theory. This energy density 
can be dubbed as dark curvature energy density, fc is the re-scaled crossover distance 
that is deflned as = r^fiR) and a prime marks differentiation with respect to the Ricci 
scalar, R. We note that in this scenario there is a n effective gra vitational constant, which is 



re-scaled by f'{R) so that G = G, 



eff 



Nozari 



2009al ). In order to compare this 



model with other alternative scenarios, it is more suitable to rewrite the normal branch of 
the Friedmann equation (16) in the following form 

Hiz) 



Ho 



(1^ 
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where by definition 



and also 



-1 



Rf"{R) + R Rf"'{R) - Hf"{R) 



'^[f{R)-Rf'{R)]-3HRf"{R) 
We note that equation of state parameter of the curvature fluid is not a constant 



(19) 



it varies actually with r edshift. To proceed 



Hu-Sawicki f{R) model ( iHu fc Sawicki 



2007 



urther, in which follo ws we consider the 



Martinelli et al. 



20091 ) given by 



f{R) = R-m 



2 



(20) 



where m^, ci, C2 and n are free positive parameters that can be 



of density parameters ( iHu fc Sawicki 



2007 



Martinelli et al. 



expressed as functions 



20091 ) ■ Here we explore the 



dependence of these parameters on density parameters defined i n our setup. To do this end 



we follow the procedure presented in Ref. ( Hu fc Sawick 



20071 ). Variation of the action 



(14) with respect to the metric yields the induced modified Einstein equations on the brane 

1 



Gap - J^-SaP - £al3 , 

where S^b, the projection of the bulk Weyl tensor on the brane is given by 



(21) 



C _(5) ^Af „ ^Rn^nS 



(22) 



and Sai3 as the quadratic energy-momentum correction into Einstein field equations is 
defined as follows 



>a/3 



4 



12 



1 

24 



(23) 
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r„/9 as the effective energy-momentum tensor localized on the brane is defined as (INozari 



2008a : 



Atazadeh et al 



2008 



Saavedra &: Vasquez 



20091) 



mz 



To^p = -mlf\R)G^p + ^yf{R)-Rf\R)\g^p + T^p + ml^^Vpf{R)-g^pUf\R)^ . (24) 

The trace of Eq. (21), which can be interpreted as the equation of motion for f'{R) , is 
obtained as 



R 



11 

6M|' 



(25) 



r, the trace of the effective energy-momentum tensor localized on the brane is expressed as 



r = m„ 



2/(i?) - Rf\R) - 3af'{R) 



(26) 



To highlight the DGP character of this generalized setup, we express the results in terms 
of the DGP crossover scale defined as = So the equation of motion for f'{R) is 



rewritten as follows 



^rl([2f{R) - Rf\R)\\9{nf'{R)y + QRf'{R)af'{R) - l2f{R)Uf\R) 



+ 



3 Ml 



11 



Rf'iR)-2f{R) + 3nf\R) 
In the next stage, we solve this equation for □/'(/?) to obtain the following solution 

'l(6i.m)-12/(i.))+|^ 



(27) 



of{R) 



± 



i(6i?/'(i?)-12/(i?) 



MfrcA 



where is defined as 



e 



2/(i?) - RfiR) 



1 



6 



llr2 



:R. 



-e, 

(28) 
(29) 



Now we introduce an effective potential 14//, which satisfies the following equation 

dVeff 



df'{R) 



(30) 
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This effective potential has an extremum at G = 

1 



2f{R)-Rf{R) 



Pr, 



re V 11 ■ 



(31) 



In the high-curvature regime, where f'{R) — 1 and ~ 1 , we recover the standard DGP 



R 



result (one can compare this result with corresponding result in Ref. (iHu fc Sawicki 
to see the differences in this extended braneworld scenario) 

2 



20071) 



R± 



6 



R 



mt 



■pr. 



(32) 



Tc V 11 

The negative and positive sign in this equation is corresponding to the DGP self-accelerating 
and normal branch respectively. In which follows, we adopt the positive sign corresponding 
to the normal branch of the scenario. To investigate the expansion history of the universe 
in our model, we restrict ourselves to those values of the model parameters that yield 



expansion histories which ar e observationa. 



function, introduced in Ref. ( IHu fc Sawicki 



y viable. We note that the Hu-Sawicki f{R) 



2007 



Martinelli et al 



20091 ). was interpreted as 



a cosmological constant in the high-curvature regime. The motivation for that interpretation 
was to obtain a ACDM behavior in the high curvature (in comparison with m^) regime. 
Here we would like to investigate f{R) models that mimic the phantom-like behavior on 
the brane in the mentioned regime. As we have pointed out previously, the phantom-like 
behavior can be realized from the dynamical screening of the brane cosmological constant. 
In this respect, we apply the same strategy to our model, so that the second term in the 
Hu-Sawichi f{R) function ( that is, second term in the right hand side of equation (20)) 
mimics the role of an effective cosmological constant on the DGP brane. Then this term 
will be screened by M- term in the late time (see the normal branch of Eq. (16)). 



In the case in which R 3> , one can approximate Eq. (20) as follows 

Ci 2 ^1 2 f ^ 



lim fiR)^ R- —m^ + ^mH — ] . (33) 
„^2/^^o C2 cl \im?J 
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During the late-time acceleration epoch, /o(-R) ^ 1 or equivalently Ro ^ and we can 
apply the above approximation. Also the curvature field is always near the minimum of the 
effective potential. So, based on Eq. (31), we have 



R+-J^R=l^p^ + 2^m2 . (34) 
Tc V 11 C2 



Since R in the f{R) function is induced Ricci scalar on the bra ne, we except crosso ver scale 



20071) they 



to affect on the constant parameters ci , C2 and m^. In Ref. ( iHu &: Sawicki] 
obtained 3m^ = R^ = that p^m is the present value of the matter density. But in our 
setup the present value of the matter density (see Eq. (32)) is given by 

2 

Tc V 11 

If we solve this equation for Rc, we find 



3m2 = Rc = l.l^r, + 6fi„ ± y 0.55fir, (o.bb^r^ + 6f]„) . (36) 

Therefore, the DGP character of this extended modified gravity scenario is addressed 
through m^. As we have argued, at the curvatures high compared with m^, the second term 
on the right hand side of equation (20) mimics the role of an effective cosmological constant 
on the brane. In this respect, the second term in the right hand side of equation (33) also 
mimics the role of a cosmological constant on the brane in the high curvature regime. With 
this motivation, we find 

£l« Z'^"-' (37) 

l.mr^ + 6fi„ ± W0.55fi^^ (o.55nr^ + 6fi^) 

There is also a relation for % as follows 

(38) 

where ^ in our setup can be calculated as follows: firstly, by using Eqs. (34) and (37), we 
find 

R+^^^ = ^p,ma-^+12Qj,, (39) 
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where p^m can be omitted through Eq. (26) to obtain 

m 18\ 



Tc V Tf. V 11^ 



(40) 



Finally, if we solve this equation for \/R , we find the following relation for ^ 



Rn 



Am 



+ 



Am 



3 1 + 



4fiA 



+ 



5aAT 



2m 



l/2\ 2 



(41) 



where m is given by Eq. (36). Note that we have set Hq and a(to) equal to unity. These 
relations tell us that the free parameters of this model are n, VLm-, ^^id /o(-R), whereas 
the latter one is constrained by Solar-System tests. In fact, experimental data show that 
f'{R) — 1 < 10~^ , when f'{R) is parameterized to be exactly 1 in the far past. To analyze 
the behavior of Wcurv , we need to specify an ansatz for the scale factor. Here we use the 



following form 



a{t) = {t^ + 



7^ 1 is a free parameter (jCai et al. 



1 - v' 



(42) 



[20071). By noting that the Ricci scalar is R = 6(f + (f )2), 
one can express the function f{R) of equation (20) in terms of the redshift z. Figure 1 
shows the variation of the effective equation of state parameter versus the redshift. As 
we see in this figure, in this class of models the curvature fluid has an effective phantom 
equation of state, Wcurv < — 1 at high redshifts and then approaches the phantom divide 
{wcurv = — 1) at a redshift that decreases by decreasing n. The main point here is that 
a modified induced gravity of the Hu-Sawicki type in DGP framework, gives a phantom 
effective equation of state parameter for all values of n. Note that all of these models reach 
asymptotically to the de Sitter phase {wcutv = —1)- As we will show via a dynamical system 
approach, this de Sitter phase is a stable phase. 



Now using equation (13) we plot the luminosity distance versus the redshift in this 
setup. The result is shown in figure 3 in comparison with other alternative models. Up to 
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the parameter values adopted here, the ADGP model has the potential to give a better fit 
with ACDM than the /(i?)-DGP model. 

As we have mentioned previously, Wcurv varies with redshift and is less than —1 which 
indicates that the curvature fluid plays the role of a phantom scalar field in the equivalent 
scalar-tensor theory As we will show in the next section, in a dynamical system approach 
the accelerating phase of the universe (g < 0) for this DGP-inspired f{R) model is stable 
if we consider the curvature fluid to be equivalent to a phantom scalar fleld ( this means 
that we set Wcurv < ^1 )) while if we set —1 < Wcurv < ~\ (equivalent to a quintessence 
scalar field ), we find that current acceleration in the mentioned universe is a transient 
phenomenon. 



4.1. The phase space of /(i?)-DGP models 

To investigate stability of the solutions presented in the previous subsection, here we 
express the cosmological equations of /(i?)-DGP scenario in the form of an autonomous 
dynamical system. For this purpose, we define the following normalized expansion variables 



In this way, equation (16) with minus sign and in a dimensionless form is written as follows 

1 + 2m = + p2 _ (^44) 

This constrain means that the phase space of this scenario in the [s - p) plane is outside of 
a circle with radius 1 which is defined as + > 1 . The autonomous system is obtained 
as follows 

/ _ 3g[g^ + (1 + 2Wcurv)p^ - 1] . 
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, _ 3p[2s2 + + s2 - 1)] 

2(s2 +J92 + 1^ 

To achieve the critical points of this system one should set s' = and p' = . Note that 
the critical points and their stability depend on the value of Wcutv ■ Here we investigate the 
stability of critical points in two different subspaces of the model parameter space where 
EoS of the curvature fluid has either a phantom or a quintessence character. In table 1, 
we see that the current accelerating phase of the universe expansion is unstable if the 
curvature fluid is considered to be a quintessence scalar field (—1 < Wcurv) and stable if it 
is considered to be a phantom field {wcurv < "l) in the equivalent scalar-tensor theory. 
This result is compatible with evolution of the deceleration parameter versus the redshift 
as has been shown in figure 4. When 2; — t- — 1, the accelerating phase (with g < 0) is stable 
for Wcurv < — 1. It is necessary to mention that whenever Wcurv = — 1, the phase space is 
ID (here the curvature fluid plays the role of a cosmolo gical constant, the same as ADGP 



model. For more details see Ref. (jChimento et al. 
trajectories of the model. 



20061 )). Figure 2 shows the phase space 



4.2. F(i?, ¥?)-DGP model 



Now we extend the previous model to an even more general case that the modified 
induced gravity is non-minimally coupled to a canonical scalar field on the brane. This 
extension allows us to see the role played by the non-minimal coupling of gravity and scalar 
degrees of freedom in the cosmological dynamics on the brane. We consider a gen eral 



coupling between grayity and scalar degrees of freedoni on the 



2008a : 



2007a : 
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2008 
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: ra & Vasauez 


2009 


Barenboim & Lvkken 


2008; 


Bamba et al. 


2OO9I) 



on the 


Drane 


Nozari 


2008b 



see for instance 



Faraoni 



2000 



Nozari 



Nozari 



In fact, inclusion of this field brings 
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the theory to reahze a smooth crossing of the phantom divide hne in a fascinating manner 



/ Tfl 

d^x^ -^F{R, 0) - ig/^- V^0V.0- 1^(0) +MlK+Cm ) , (47) 



The action of this general model is given as follows f lNozarill2009bl ) 

V 2 

where the first term shows the usual Einstein-Hilbert action in the 5D bulk. The second 
term on the right hand side is a generalization of the Einstein-Hilbert action induced on the 
brane. This is an extension of the scalar-tensor theories in one side and a generalization of 
/(i?)-gravity on the other side. K is the trace of the mean extrinsic curvature on the brane 
in the higher dimensional bulk, corresponding to the York- Gibbons- Hawking boundary 
term. We call this model as F{R, (/))-D GP scenario. N ote that from the above action, the 



crossover scale takes the following form (INozari 



2009ah 



2M| 



F'{R,ip)=r, F'{R,ip) 



(4J 



where as usual 



and a prime denotes a differentiation with respect to R. Since 



PGP s cenario accounts for embedding of the FRW cosmology at any distance scale (INozari 



2007bl ). we start with the following line-element 



ds'^ = q^ydx^dx" + t)dif' = —n^{y, t)dt^ + a^{y, t)'jijdx^dx^ + b'^{y, t)dy'^ 



(49) 



where is a maximally symmetric 3-dimensional metric defined as = 5ij + k^^^. 
Also k = —1,0, 1 parameterizes the spatial curvature and = XjX*. We choose the gauge 



^ We note that the normal branch of the pure DGP scenario ( which has the capability 
to describe the phantom like effect), cannot realize crossing of the phantom divide line 
without introducing a quintessence scalar field (a canonical field) on the brane. Introduction 
of a quintessence field on the brane, brings the theory to realize this interesting feature 



( IChimento et al 



20061 ) . For F{R, ip)-DGP scenario it i s natural to expect re alization of the 



this feature due to wider parameter space in this case (INozari 



2008a 



2009ah . 
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b'^{y,t) = 1 in the normal Gaussian coordinates. The cosmological dynam ics on the bran e 



in this model can be described by the fo 



Atazadeh et al 



2008 



Saavedra &: Vasquez 



owing Friedmann e quation (see 



Nozari 



2008a 



2009 



Nozari 



2009ah 



3mlF'{R, ip) 



Jot 



Jot 



mlF'{R,ip) 



6M| 



(50) 
and 



where e = ±1 shows two different embedding of the brane in the bulk, g 
^0 = 3^^^ — fl'^ + k^ci^ is a constant with respect to y ( with a' = see ( Dick 2001 ) for 
more detailed discussion on the constancy of this quantity. Total energy density and pressure 

Pm + P0 + p^'^"^'^) + Pa and = Pm + P<f> + p^^^^""^ + pA respectively. 



are defined as p 



(tot) 



The ordinary matter on the brane has a perfect fluid form with energy density pm and 
pressure Pm, while the energy density and pressure corresponding to non- minimally coupled 
quintessence scalar field and also those related to curvature fluid are given respectively as 
follows 



1 . rIF 

-02 + nV(0) - 6^i/0 



(51) 



2n^ d(f) 



^0) + 4^^^0+ 2 d^F^^- 
n d(f) V? V? dcjP' 



y=0 



(52) 



P 



(curv) 



P 



(curv) 



2/1 



F{R,(f)) - RF'{R,, 



m;i2RHF"{R, <f)) + RF"{R, 0) + R'F"'{R 



3RHF"{R 
1 



Also by definition p\ = Ricci scalar on the brane is given by 



R 



k 



1 



6- 



, \ 2 ah 
6(- -6 — 

.a/ an 



F{R, (f)) - RF'{R, , 



(53) 
• (54) 



We set n{0,t) = 1 where y = is chosen to be the location of the brane. By neglecting the 
dark radiation term in equation (23), we find 
SvrG 



3-(p.+P.V-™^) + | + ^ + eJ^ + ^ 



SnG , , A" 

^(p™ + P^ + p(^"™))+ 3 , (55) 
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where G = Geff = ^^^2p,(^ji ■ By adopting the negative sign, we have 



H 



Comparing this equation with the following Friedmann equation 



(56) 



SttG 



(57) 



we are led to conclude that the screening effect on the cosmological constant is modified by 
F'{R, if) as follows 

^^^^ ~ 3 r,F'{RM ■ ^^^^ 
To see how this model works, we consider an explicit form of F[R^ if) as F{R, if) = 
^(1 — E(t>'^)[R — (1 — n)C^(R/C^)^] where C is a suitably chosen paranaeter ( s ee for instan ce 



Capozziello 



2003 



Sotiriou 



2010 



Noiiri fc Odintsqv 



2003 



2004 . 



2005 



2006 



Noiiri et al. 



2007 



Noiiri fc Odintsov 



Starobinsky 



2008 



Bamba et al. 



20071 ) and (ICapozziello et al. 



2008 



2008 



Carroll et al 



2004 



Amendola et al. 



20071 : 



20101) • For a spatially flat FRW geometry. 



the Ricci scalar is given by i? = 6^ + 6(^)^. One can investigate variation of the effective 
dark energy density with respect to parameters ^, t and n to see t he status of phantom 



2009al ). this model accounts 



mimicry in this model. As has been shown in Ref. (iNozaril 
for, and modifies the phantom like behavior on the brane. Using the conservation equation 
Peff + 3i7(l + u}p,ff)pp,ff = 0, fulfillment of the condition Wefj < — 1 requires the constraint 



H F'{R,4>) 



2009a|). In this situation, this model has the potential to realize the 



phantom-like behavior and smooth crossing of the phantom divide line by the effective 
equation of state parameter. Now, in order to compare the F{R, if)-DGP model with 
ACDM and also with other alternative models, we study expansion histories of these models 
based on the variation of their luminosity distances versus the redshift. The evolution of 
the cosmic expansion in the normal branch of this F(R,(y9)-DGP setup is given by 

= fi„(l + Z)^ + + ;2)3{1+--) + fi,„,,(l + ^)3{l+-c„..) + + 2fi,, 

-fin 
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- 2 yn., ^Jftm{l + + Z)3(l+-^) + + ^)3(l+c.™..) + Q^^ . (59) 

Note that by definition, = ^t'^h^f'^r ^) ^^cl we use the normalization F'(_R, y9)|(2=o) — 1- 
This model has the disadvantage that it contains more free parameters than ACDM and 
even ADGP. However, this larger parameter space provides, at least theoretically, more 
facilities to realize some interesting features. The parameter space of this model can be 
constraint at 2; = so that 



(60) 



Since VLr^ fi^ + VL^p + Vtcurv + + should be p ositive, there are two possibiliti es 



for constrains on the parameters of t 



2004 



Shtanov 



2000 



lis m odel as follows (jSahni fc Shtanov 



Lue &: Starkman 



2003 



Sahni 



2OO4I ) 



(61) 



+ fiA + 2Vl]. = 1. 



(62) 



In our forthcoming analysis, we consider the first one of these constraints. Finally, using 
equation (13), we compare the luminosity distance - redshift relation of this model with 
other proposed scenarios. The results are shown in figure 3. In comparison with ADGP and 
/(-R)-DGP, the F{R, <^)-DGP scenario has the potential to have a better fit with ACDM. 
The parameters adopted in this model are Vt^ = 0.27, = 0-01 5 = 0.1 , Qcurv = 0.13 
and ^A = 0.7. 



4.3. ^(6", <^)-DGP model 



In the F{G, (y9)-DGP setup, the curvature corrections are taken into account via 
incorporation of the Gauss-Bonnet invariant term in the brane part of the action. The 
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Gauss-Bonnet invariant should be considered essentially in the bulk action. But, if there is a 
coupling between scalar degrees of freedom on the brane and the mentioned invariant term, 
then it is necessary to consider its contribution in the field equations on the brane too. This 
is the main reason for incorporati ng scalar field in the F(G,(fi) setup. We start with th e 



action of this scenario as follows (iNozari et al 



2009 



Nojiri et al. 



2007 



Bamba et al. 



20101) 



M3 



mz 



^-^R + MlK + FiG,^) + Cr, 



The Scalar- Gauss-Bonnet term in this action is defined as (INozari et al. 

1 



20091 ) 



(63) 



By definition, the Gauss-Bonnet invariant G{R) is given by 

G{R) = R^- 4R^.R^- + R^^^pR^"''^. 
In order to discuss about cosmological dy namics on the bra ne, one can achieve the 



(64) 



generalized Friedmann equation as follows (INozari et al 



20091 ) 



3m2 



+ 



2r2 



+ e. 



1 fp^ + pG^ 



4^4 



3m2 



The energy density corresponding to the Gauss-Bonnet term is defined as 



(65) 



P 



02 + ^(0) _ 2A^f'{cP)H' 



(66) 



and the corresponding pressure is defined as 



p^-' = _ v(0) + 8-(//V) + 160/'(0)/^^ 



(67) 



where a dot marks differentiation with respect to the cosmic time. One can express the 
Friedmann equation (65) for the normal branch of this DGP-inspired scenario and for an 
spatially flat brane in a dimensionless form as follows 
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- 2^^^nm{l + Zf + Qgs(1 + Zf^^+'^OB) + (68) 

where by definition Q,gb — f^Po'^^^- ^^^^ red-shift 2; = 0, equation (68) can be expressed 
as 

1 = (yaT+TwTa^ - \ (69) 

By choosing the positive sign of the square root, this constrain equation yields 

n^ + ^GB-'^V^^^l. (70) 

Note that the general relativistic limit can be recovered if we set Q^c — (or M5 = 0). In 
this case equation (70) implies Jl^ + VIgb — 1- Finally, the luminosity distance-redshift of 
this model is compared with other proposed alternatives in figure 3. 



5. Some other probes 

To have more complete comparison between the proposed models, in this section we 
study some other features of these models. The current universe is accelerating and this 
can be realized by investigating the deceleration parameter. The deceleration parameter 
q — —-^a can be expressed as 

H[z) 

where a prime marks differentiation with respect to z. We study variation of q versus 
the redshift in each of the previously proposed models. For this purpose, we use H[z) 
as given by equations (9), (12), (18), (59) and (68). Figure 4 compares the deceleration 
parameter in each of the proposed scenarios. All of these scenarios explain the late-time 
cosmic speed-up in their normal DGP branches, but the redshift at which transition to the 
accelerated phase occurs are different in these scenarios. The ADGP model transits to the 
accelerated phase at 2; ~ 0.8 much sooner than other scenarios. The F{R, ip)-DGP model 
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transits to the accelerated phase at z ~ 0.4 much later than other alternatives. We note 
however that the exact value of the transit redshift depends on the choice of the parameters 
in each model. Within the parameters values adopted here, the result are shown in figure 4. 
As another important result, we note that late-time acceleration in the f{R)-DGP universe 
is a stable phenomenon since we considered the curvature fluid to be a phantom scalar 
field. Also the acceleration phase is a transient phenomenon for F{G, 99)-DGP model if 
we set > —1. Our inspection shows that cosmological dynamics in the normal branch 
of /(i?)-DGP scenario with a quintessence curvature fluid is very similar to cosmological 
dynamics in the normal branch of F{G, 9?)-DGP scenario for wgb > — 1- 

The next probe is the age of the universe in each of these scenarios. The age of the 
universe at a given cosmological red-shift is given as follow 



where H{z) is given in each of the proposed scenarios. In figure 5 the age of the universe 
is compared in five alternative scenarios. Within the parameter values adopted here, the 
age of the universe in /(it!)-DGP is larger than ACDM, but this age in F{G, ip)-DGP is 
smaller than ACDM one. 



In this paper, observational status of the ACDM , ADGP and some extended 
phantom-like cosmologies as alternatives for dark energy are studied. We compared these 
models with each other by comparing variation of their luminosity distances versus the 
redshift and also the age of the universe in each of these scenarios. The common feature 
of the mentioned DGP-inspired scenarios is possibility of realization of the phantom-like 
behavior in their normal DGP branches. In the extended scenarios, we focused on the 




(72) 



6. Summary 
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role played by new ingredients such as the curvature effect, modification of the induced 
gravity on the brane and even the existence of a quintessence field non-minimally coupled 
to modified induced gravity on the brane. Since all of our analysis are preformed on the 
normal DGP branch of each scenario, there is no ghost instabihties in these setups. As an 
example and to have an intuition about stability of the solutions, we analyzed stability of 
cosmological phases of /(i?)-DGP scenario within a dynamical system approach. We have 
analyzed the phase space of the normal branch of this model where curvature fluid plays 
the role of a phantom or quintessence scalar field in the equivalent scalar-tensor theory. We 
have shown that the matter dominated phase of this model is a repeller, unstable point. 
However, there is a de Sitter phase which is an attractor point if curvature fluid plays the 
role of a phantom scalar field in the equivalent scalar-tensor theory. For a quintessence 
field plying the role of the curvature fluid, the de Sitter phase is an unstable, saddle point. 
Up to the parameters values adopted in this paper, F{R, (/9)-DGP scenario mimics the 
behavior of the ACDM scenario more than other proposed models. As another important 
outcome, all of these scenarios explain the late-time cosmic speed-up in their normal DGP 
branches, but the redshift at which transition to the accelerating phase occurs are different 
in these scenarios. For instance, the ADGP model transits to the accelerating phase at 
z ~ 0.8, much sooner than other scenarios. However, the F{R, (p)-DGP model transits to 
the accelerated phase at z ~ 0.4, much later than other alternatives. Nevertheless, the 
exact value of the transit redshift depends on the choice of the parameters in each model. 
Finally, we compared the age of the universe in each of the proposed scenarios. Within the 
parameter values adopted in this paper, the age of the universe in the /(it!)-DGP model is 
larger than AGDM, but this age in F{G, (p)-BGP is smaller than ACDM one. 
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Fig. 1. — Wcurv versus the redshift for a DGP-inspired f{R) model. The f{R) is chosen to 
be the Hu-Sawicki model. The selected free parameters in plotting this figure are chosen to be 
= 0.27 = 0.93 jQj-c = 0.01 and fQ{R) — l = 10^^ . As this figure shows, in this class of models 
the curvature fluid has an effective phantom equation of state with Wcurv < — 1 at high redshifts. 
This effective equation of state parameter approaches the phantom divide line {wcurv = — 1) at a 
redshift that decreases by decreasing the value of n. 



- 33 - 




I / 

/ / ^-"-v \ \ \ 
/ Z^^-'-.x \ \ \ 

/ ^^^v \ \ \ \ 

/ ^♦-'^ V \ \ \ \ 
^^^^ \ \ \ \ \ 

^^■^ v\ \ \ \ \ t 

\ \ \ 

W\\ \ \ \ \ t 

\ WV\ \ \ \ t 
W'^^ V\ \ \ \ 
\ 

1 T" 

1 2 
S 




t t t t 



^ N V \ \ \ 1 ■ ■ 



Fig. 2. — The phase space of the normal branch of DGP inspired f{R) models in which curvature 
fluid plays the role of a phantom scalar field (left) and the quintessence scalar field (right) in the 
equivalent scalar-tensor theories. Point A in these two cases is the matter dominated phase and it 
is a repeller, unstable point. Whereas, point B that represents the de Sitter phase, is an attractor 
point in the first case and a saddle point in the second case. 
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Fig. 3. — Luminosity distance versus the red-shift for five alternative scenarios proposed in this 
paper. In this figure we have set {Qm, ^r-c) = (0.27, 10~^) for all of the models and {flip, J^a, ^curv) = 
(0.1, 0.7, 0.13) for the F{R, ip)-DGP setup, = 0.93 for the ADGP, Qcurv = 0.93 for the /(i?)-DGP 
setup and flcB = 0.93 . We calculated these quantities using the constraint equations attributed 
to each model. Up to the parameters values adopted here, F{R, (/?)-DGP scenario is closer to the 
ACDM than other proposed models. We note however that this result depends strongly on the 
adapted parameters values. 
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Fig. 4. — The deceleration parameter q{z) versus the red-shift for five alternative scenarios dis- 
cussed in this paper. The model parameters are the same as adopted in figure 3 and Wcuw here is 
adopted quintessence like in the left and phantom like in the right. Note that current acceleration 
in the /(i?)-DGP universe in which curvature fluid play role of a quintessence is a transient and 
unstable phenomenon but if it act as a phantom field, it result in a stable accelerated phase for the 
universe. 
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Fig. 5. — The age of the universe (in units of the inverse Hubble parameter) versus the red- 
shift for five alternative scenarios discussed in this paper. The model parameters are the same as 
adopted in figure 3. With proposed parameter values, the /(i?)-DGP model has the largest and 
the F{G, (p)-DGF has the smallest age. 
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Table 1: Eigenvalues and the stability properties of the critical points. 

points {s,p) eigenvalues for Wcurv < ~1 for (~1 < 'Wcurv < 0) 

A (1,0) ^- ^wcurv ^ 1^ unstable unstable 

B (0,1) (3^^ 3(i+^c...) -) g^j^^^g unstable 



